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This letter reports a measurement of the cross section of prompt isolated photon pair production in 
pp collisions at a total energy yfs = 1.96 TeV using data of 5.36 fb _1 integrated luminosity collected 
with the CDF II detector at the Fermilab Tevatron. The measured cross section, differential in basic 
kinematic variables, is compared with three perturbative QCD predictions, a leading order (LO) 
parton shower calculation and two next-to-leading order (NLO) calculations. The NLO calculations 
reproduce most aspects of the data. By including photon radiation from quarks before and after 
hard scattering, the parton shower prediction becomes competitive with the NLO predictions. 
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new phenomena, such as new heavy resonances [2j , extra 
spatial dimensions 0, 0] or cascade decays of heavy new 
particles [j| . Precise measurements of the diphoton pro- 
duction differential cross sections for various kinematic 
variables and their theoretical understanding are thus 
very important for these searches. Diphoton production 
is also used to check the validity of perturbative quantum 
chromodynamics (pQCD) and soft-gluon resummation 
methods implemented in theoretical calculations. Dipho- 
tons are expected to be dominantly produced by quark- 
antiquark annihilation qq — > 77 and also in kinematic 
regions with high gluon luminosity, especially at low in- 
variant mass, by gluon-gluon fusion gg — » 77 through 
a quark loop diagram. Prompt photons may also result 
from quark fragmentations in hard scattering, although a 
strict photon isolation requirement significantly reduces 
the fragmentation contributions. 

Diphoton measurements have been previously con- 
ducted at fixed-target @ and collider experiments @- 
[Toj . The most recent measurements [^, [13] were com- 
pared with the same pQCD calculations examined in the 
present work and large discrepancies were found between 
the data and a LO matrix element calculation supple- 
mented with a parton shower model, suitable for simula- 
tion of the backgrounds in searches of a low mass Higgs 
boson and of new phenomena. This work shows that the 
inclusion of photons radiated from initial and final state 
quarks drastically improves the comparison of the parton 
shower calculation with the data. 

The reported measurement was conducted using data 
of total integrated luminosity 5.36 fb _1 collected with 
the Collider Detector at Fermilab (CDF) [ll| at the 
Tevatron pp collider. CDF is composed of a cen- 
tral spectrometer inside a 1.4 T magnetic field, sur- 
rounded by electromagnetic and hadronic calorimeters 
and muon detection chambers. The inner spectrometer 
measures charged particle tracks with a transverse mo- 
mentum (pt) precision of Apx/Px = 0.07%(GeV/c)~ 1 . 
The central calorimeters cover the region |?/| < 1.1, 
with an electromagnetic (hadronic) energy resolution of 
a(E T )/E T = 13.5%/ ^/E T (GeV) © 1.5% (<j(E t )/E t = 
50%/i/i?T(GeV) © 3%) and a tower segmentation of 
A77 x Acf> ~ 0.1 x 15°. Photons are reconstructed in clus- 



ters of up to three towers [12j. \ 2 criteria are imposed 



on the profile of the shower to match expected patterns. 
Two main cuts are applied: (i) the photon transverse en- 
ergy is required to be Et > 17 GeV for the first photon 
in the event and Et > 15 GeV for the second photon; 
(ii) the calorimeter isolation energy in the isolation cone 
around each photon [l3| is required to be less than 2 
GeV. 

The background from 7+jet and dijet events, where 
one or two jets are faking a photon, is subtracted with a 
method using the track isolation as the discriminant be- 
tween signal and background [3] ■ It is based on the sub- 
stantial difference of the track isolation distribution for 



signal photons (nearly exponential) and for background 
photons (nearly flat). The advantages of this method are 
that (i) it has little sensitivity to multiple interactions 
in the colliding beams, so that the signal-background 
separation does not degrade at high instantaneous lu- 
minosity, and (ii) it has high efficiency and good track 
momentum resolution, implying minimal degradation of 
the signal-background separation due to instrumental ef- 
fects. The signal fraction is determined by summing the 
probabilities of an event to be pure signal, pure back- 
ground, or mixed photon pair. These probabilities are 
obtained by solving a 4 x 4 matrix equation using the 
observation value (0 or 1) for all four combinations of 
the leading or sub-leading photon having track isolation 
below or above 1 GeV/c as an input. The matrix is 
constructed from the -&r-dependent efficiencies of signal 
and background photons passing the track isolation cut. 
A threshold cut of 1 GeV/c is determined by maximizing 
the separation between signal and background. The ef- 
ficiencies are determined from Monte Carlo (MC) 7+jet 
and dijet samples, which are produced using the PYTHIA 
event generator [l5j. PYTHIA events are fully simulated 
through the detector and trigger and are reconstructed 
with the CDF II simulation and reconstruction software 
[3] . With this matrix technique the full correlations be- 
tween the two photons in the event are properly taken 
into account. Tests are made for underlying event con- 
tributions in complementary cones to the photon recon- 
struction cone [17j | and also using isolated tracks in dijet 
events. The systematic uncertainty in the signal fraction 
with this method is of the order of 15-20%. 

The diphoton production cross section differential in a 
kinematic variable is obtained from the histogram of the 
estimated signal in the selected variable. The average 
cross section in a bin of the variable is determined by di- 
viding the bin content by the trigger efficiency, the dipho- 
ton selection efficiency and acceptance, the integrated lu- 
minosity and the bin size. The diphoton trigger efficiency 
is derived from data [lj. It is consistent with 100% over 
all of the kinematic range with a flat uncertainty of 3%. 
The selection efficiency is determined from data and MC 
with an iterative method. In the first pass the efficiency 
is determined from a fully simulated and reconstructed 
PYTHIA diphoton MC sample by dividing the number of 
events passing all selection cuts by the number of events 
passing only the kinematic cuts on the photon E^, rj, 
angular separation and isolation at the event generation 
level. The efficiency denominator is corrected for the "un- 
derlying event" from collision remnants which make the 
efficiency obtained from pythia too high by removing 
events from the denominator through the isolation cut. 
This correction is derived by running pythia with and 
without underlying event and amounts to a constant fac- 
tor of 0.88 per event. A flat 6% uncertainty in the selec- 
tion efficiency (3% per photon) accounts for possible in- 
accuracies in the pythia model for the underlying event. 
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The signal events of the data are corrected for the prelim- 
inary efficiency. The data are then used to reweight the 
pythia events and obtain a more accurate representation 
of the true diphoton distribution. The efficiency is deter- 
mined using the reweighted PYTHIA sample and corrected 
for luminosity dependence, derived from a comparison of 
the vertex multiplicity distribution in data and pythia 
MC Z° — > e + e~ events. The systematic uncertainty in 
the efficiency resulting from the luminosity dependent 
correction grows linearly from 1.8% for Et < 40 GeV 
to 3% at Et = 80 GeV and remains constant above this 
point. Finally, a 6% constant uncertainty comes from the 
Tevatron integrated luminosity (l8j . 

The Z° — > e + e~ sample is used for calibration by ap- 
plying a "diphoton-like" event selection, i.e. by impos- 
ing a diphoton selection with the same trigger but al- 
lowing for a track associated with each of the two elec- 
tromagnetic objects in the event. The electromagnetic 
energy scale in data and MC is corrected by tuning the 
Z° —> e + e~ mass peak to the world average and a 
systematic uncertainty from this correction is estimated 
to grow linearly from at Et < 40 GeV up to 1.5% 
at Et — 80 GeV and remain constant above this point. 
The difference in the photon identification efficiency be- 
tween data and MC is estimated from the Z° — > e + e~ 
sample [l[ and added as a systematic uncertainty to the 
measurement. All systematic uncertainties in the cross 
section measurement are added in quadrature. 

The results of this measurement are compared with 
three theoretical calculations: (i) the fixed NLO predic- 
tions of the DIPHOX program [20| including parton frag- 
mentations into photons [2l| , (ii) the predictions of the 
RESBOS program [22[ where the cross section is accurate 
to NLO, but also has an analytical initial state soft gluon 
resummation, and (iii) the predictions of the pythia 
program (l5| which features a realistic representation of 
the physics events by including parton showering, Initial 
(ISR) and Final State Radiation (FSR) and an underly- 
ing event model. Diphoton events were selected from an 
inclusive 7+X pythia sample (X=7 or jet), thus includ- 
ing the qq — > 77 and gg — > 77 processes (56%) as well 
as the qq ->■ 577 ISR , gq -> ?77i S r and 99 -> 977fsr pro- 
cesses (44%). This type of calculation effectively resums 
the cross section for gluon and photon radiation both in 
the initial and the final state. All calculations are subject 
to the experimental kinematic and isolation cuts, diphox 
accounts for the gg — » 77 process in LO only. The predic- 
tions of RESBOS are restricted to the invariant mass range 
from 2m,b — 9 GeV/c 2 to 2m t = 350 GeV/c 2 , where nib 
and Tfit are the masses of the bottom and top quarks, re- 
spectively. NLO theoretical uncertainties are estimated 
by varying the fragmentation (in diphox only), renor- 
malization and factorization scales up and down by a 
factor of two relative to the default scale /i = M/2 of 
diphox and (j, = M of RESBOS, and for the NLO PDF 
uncertainties (in both DIPHOX and RESBOS) by using the 



20 CTEQ6.1M eigenvectors [H]. 

Figure [T] shows the comparison between the measured 
and predicted diphoton distributions: the diphoton in- 
variant mass M, the diphoton transverse momentum Pt 
and the difference Acf> between the azimuthal angles of 
the two photons in the event. While the pythia direct 
calculation (77) fails to describe both the scale and shape 
of the data, including radiation brings the prediction in 
fair agreement with the data. In particular, radiation 
makes the Pt and A<fi distributions harder because of 
the presence of at least one hard jet in the final state 
of events in which one photon originates from radiation. 
The mass distributions show a reasonable agreement with 
the data for all predictions above the peak at 30 GeV/c 2 , 
particularly in the region 80 GeV/c 2 < M < 150 GeV/c 2 
relevant to searches for the Higgs boson [![. However, all 
predictions underestimate the data around and below the 
peak. In the Pt spectrum all predictions underestimate 
the data in the region between 20 and 50 GeV/c, a fea- 
ture also observed in the earlier measurements [9j, |20( • For 
Pr<20 GeV/c, where soft gluon resummation is most im- 
portant, only the RESBOS prediction describes the data. 
Discrepancies between data and theory are most promi- 
nent in the comparison of the measured and predicted 
distributions of A<f>. In this case all three predictions 
fail to describe the data across the whole spectrum. Ap- 
proaching A(j) = 7r, where soft gluon processes are ex- 
pected to manifest, the RESBOS prediction agrees better 
with the data. In the range 1.4 rad<Ac/(><2.2 rad only the 
pythia prediction describes the data and remains closest 
to the data down to 1 rad. In the low A</> tail, which cor- 
responds to the region of low M (<50 GeV/c 2 ), all three 
predictions are lower than the data, although the diphox 
prediction, by explictly including non-perturbative frag- 
mentation, lies closer to the data for A</><1 rad. 

In summary, the diphoton production cross section, 
differential in kinematic variables sensitive to the reac- 
tion mechanism, is measured using data corresponding 
to an integrated luminosity of 5.36 fb _1 collected with 
the CDF II detector. The high statistics of the measured 
sample allows for a higher precision scan over a much 
more extended phase space than previous measurements. 
The overall systematic uncertainty is limited to about 
30%. The results of the measurement are compared with 
three state-of-the-art calculations, applying complemen- 
tary techniques in describing the reaction. All three cal- 
culations, within their known limitations, reproduce the 
main features of the data, but none of them describes 
all aspects of the data. The inclusion of photon radia- 
tion in the initial and final states significantly improves 
the PYTHIA parton shower calculation (see the left-hand 
panels in Figure [1]), which is suitable for background sim- 
ulations in searches for a low mass Higgs boson and new 
phenomena. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
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FIG. 1: The measured differential cross sections compared with three theoretical predictions discussed in the text. The left 
windows show the absolute comparisons and the right windows show the fractional deviations of the data from the theoretical 
predictions. Fractional deviations for pythia refer to the 77 + 7j calculation. Note that the vertical axis scales differ between 
fractional deviation plots. The comparisons are made as functions of the diphoton mass (top), transverse momentum (middle) 
and azimuthal angle difference (bottom). The shaded area around the data points indicates the total systematic uncertainty 
of the measurement. 
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